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An analysis of the sensitivity of gaussian and mexican hat wavelet family filters to the 
detection of point sources of ultra-high energy cosmic rays was performed. A source 
embedded in a background was simulated and the number of events and amplitude of 
this source was varied aiming to check the sensitivity of the method to detect faint 
sources with low statistic of events. 
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1. Introduction 

The origin of the Ultra-High Energy Cosmic Rays (UHECRs) is still an unsolved 
problem in astroparticle physics. A model for acceleration of cosmic rays, firstly 
devised by Hillas^ predicts that cosmic rays are accelerated to the highest energies 
(above 10 18 eV) by electromagnetic fields of astrophysical objects. Therefore, the 
identication of possible astrophysical sources of UHECRs is possible by analysing 
the arrival directions of the cosmic rays. The correlation of the positions of point- 
like astrophysical objects (point sources) with the arrival directions of cosmic rays 
defines a small scale anisotropy. In the search of point-like sources it is a common 
procedure to convolve the sky maps containing arrival directions of cosmic rays 
with mathematical functions (the kernel of the convolution operation) aiming to 
optimize the signal to noise ratio. In this work it is studied the performance of some 
kernels of the Mexican Hat Wavelet Family (MHWF) to identify point sources of 
cosmic rays, and disentangle genuine signals from the background. 

2. Wavelets 

Wavelets are defined as mathematical functions belonging to the L 2 space. They can 
be thought as localized wave-like oscillating functions which can be operated with 
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a given signal and provide information about it. The continuous wavelet transform 
(CWT) in two dimensions may be formally written as 

$(s,r 1 ,r 2 )= J f f(t,u)*l TuT2 (t)dtdu, (1) 

where s (s > 0, s g R) is the scaling factor and n and T2 (t.; G M.) are the 
translation parameters. So, the CWT decomposes a function f(t,u) in a basis of 
wavelet fy StTltT2 (t,u). One can scale and translate a "mother-wavelet" ^ and obtain 
a wavelet \l/ SiTl . T2 (t, u), as follows: 

1 ft — T-l U — To \ 
* S , T1 , T3 (*,U) = -^W — _J:,__»J. (2) 

The Mexican Hat Wavelet Family (MHWF), introduced by Gonzalez-Nuevo et 
al®, and its extension on the sphere have been wid ely us ed to detect point sources 
in maps of cosmic microwave background radiationSEEl, due to the amplification 
of the signal-to-noise ratio (SNR) when transiting from real to wavelet space. The 
MHWF is obtained by successive application of the laplacian operator to the two- 
dimensional gaussian 4>(x). A generic member of this family, of order n, is: 

(— 1 \ n 

*n(*0 = ~^^ 2n m- (3) 

3. Celestial Maps 

Celestial maps are pixelations of the celestial sphere taking into account the angular 
resolution of the experiment. The events map is a celestial map representing the 
arrival directions of cosmic rays in a suitable coordinate system. Due to intrinsic 
limitations of detector, every event detected is convolved with a probability related 
to the angular resolution of the detector, which means that there is a point spreading 
function (PSF) associated to the detector. 

The convolution of celestial maps with filters is given by: 

Y J ,M(m(r k ,f j ) 

where M(j) is the number of cosmic rays within the pixel of index j, in the direction 
fj. $(r, ro) is the used filter and rl is the position vector representing the point 
where the integral is being calculated. 



4. Analysis Procedure 

This work is an extension of previous oned@J2 The simulated detector has two sites, 
one located in the southern hemisphere (36° S and 65° W), and the other in the 
northern hemisphere (38° N and 102° W), seven times larger than the one in the 
south, implying on a flux of cosmic rays seven times greater. The acceptance law has 
the form sin8cos8, where 0° < 9 < 60° is the zenith angle. It was also considered 
the case of an ideal detector with uniform exposure and full sky coverage. 
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It was simulated a point source located at (l,b)=(320°,30°) (galactic coordi 
nates). This source was modeled by a gaussian: 

A fx 2 



ex P I ) ' ^ 



and was embedded in a background consisting on a superposition of four differ- 
ent patterns of arrival directions of cosmic rays. The simulated patterns were: 
(i) an isotropic distribution of events; (ii),(iii) dipoles, modeled according to 

Q(u) — ^1 + aD.uj , where $ is the flux of cosmic rays, $o is related to the 
isotropic flux, < a < 1 is the amplitude of the dipole (a = 0.07 for (ii) and 
a = 0.005 for (hi)), D is the vector which points towards the dipole ((l,b) — 
(0°, 0°)) for (ii) and (I, M = (166.5°, -29°)) for (hi)) and u is a unit vector pointing 
in an arbitrary direction^ (iv) several sources modeled according to equation [5l in 
the directions (l,b): (0°,0°) [a = 7.0°, A = 1.00], (320°, 90°) [a = 1.5°, A = 0.05], 
(320°, -40°) [a = 0.5°, A = 0.01], (220°, 10°) [a = 3.0°, A = 0.05], (100°, -70°) 
[cr = 2°, A = 0.10], (240°, 50°) [a = 20°, A = 0.05], (350°, -80°) [a = 6.0°, 
A = 0.005], (100°, 50°) [cr = 30°, A = 0.50], (140°, -40°) [a = 4.0°, A = 2.00] and 
(60°, 50°) [cr = 3.0°, A = 0.02]. This last background pattern was included in the 
simulation because unknown sources might be present during an analysis tagged 
on a given source, and their effects must be evaluated. From this combination of 
source (signal) and background patterns (noise), it was obtained the events map. 

The events map resulting from the sum of the background patterns and the 
simulated source was convolved with gaussian and MHWF (orders 1, 2 and 3) 
filters. The amplification of the SNR (A) is calculated by the expression 

A = ^, (6) 
w /<jo 

where wq is the value of the central pixel associated to the source in the non-filtered 
source map, Wf is the corresponding value in the filtered source map, ctq is the root 
mean square (RMS) of the non-filtered background map and 07 is the RMS of the 
filtered background map. 

Aiming to study the impact of the number of events from the source and its 
intensity upon the filter, different number of events (N ev t) were simulated in the 
direction of the source, ranging from 10 events up to 1000. The amplitude of the 
gaussian source (A) was also varied, from 10~ 4 to 1. 



5. Results 

In figures [T] and [2] it is shown the maximum amplification of the SNR, as a function 
of the corresponding scale. In the case of figure [1] N evt is fixed and A is varied. For 
figure [21 A is fixed and N evt from the source is varied. 

From figure [1] it can be noted that the uniform exposure acts like a constraint 
for A, and that the gaussian has a slightly better performance than the MHWF 
filters. The maximum amplification for nonuniform exposure is clearly achieved by 
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Fig. 1. Maximum amplification of the SNR (A) as a function of the corresponding scale (<Tf). The 
graphs displayed refer to a fraction of events between the source and the background of 6 X 10~ 6 
(top) and 6.25 X lO" 4 (bottom). 
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Fig. 2. Maximum amplification of the SNR (A) as a function of the corresponding scale (07). The 
graphs displayed refer to an amplitude A=10 -4 (top) and A=l (bottom). 
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using MHWF niters, whereas the gaussian filter has amplification close to 1, which 
means no amplification. Comparing the two graphs in figure [1] in the bottom graph 
the amplification is greater, which seems reasonable since the number of events in 
this case is 100 times greater. 

In figure [2] it can be seen the behavior of the filters when A is varied. For the 
uniform exposure the amplification of both the gaussian and the MHWF filters are 
low, but the gaussian has a slightly better performance. Comparing the top and 
bottom graphs in figure [21 it is clear that the amplification achieved by the filters 
is proportional to N evt . Also, when there is an acceptance, the gaussian filter does 
not provide a good amplification of the SNR, which can be achieved by using the 
MHWF. 

6. Conclusions 

In this work it was analyzed the performance of the gaussian and the MHWF 
filters to detect point sources of cosmic rays embedded in a non uniform back- 
ground, whose features are modulated both by the acceptance of the detector and 
the background patterns imposed to the incoming particles. Some parameters from 
the source such as the amplitude A and the number of events N evt were varied, and 
the effects of theses changes on amplification of the SNR was studied. 

The trivial conclusion is that the amplification achieved by these kernels is 
proportional to the source intensity parameters (A and N evt ). It is interesting to 
notice that for a realistic case, i. e., a source with low amplitude and only a few 
events coming from its direction, the amplifications achieved are low for both filters. 
The MHWF filters are more robust to these parameters and can provide a greater 
amplification of the SNR even if N evt is small. 

Regarding the contribution of the acceptance of the experiment for the detection, 
it can be clearly seen that for uniform exposure the amplifications are smaller. In 
this case, the gaussian filter provides a slightly better amplification compared to 
the MHWF. However, for a more realistic case taking into account the nonuniform 
exposure of the detector, the MHWF filters always achieve a greater amplification. 
Also, they are more robust to low statistic of events, which makes them particularly 
useful for cosmic ray studies. 
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